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We study the semileptonic decay Bs —>■ which is induced by b —>■ ui~vi transitions at 

the quark level. We take into account the standard model (SM) operator from W-boson exchange 
as well as possible extensions from physics beyond the SM. The secondary decay —>■ Kn can 

be used to study a number of angular observables, which are worked out in terms of short-distance 
Wilson coefficients and hadronic form factors. Our analysis allows for an independent extraction of 
the Cabibbo-Kobayashi-Maskawa matrix element |Kii,| and for the determination of certain ratios 
of Bs —>■ K* form factors. Moreover, a future precision measurement of the forward-backward 
asymmetry in the Bs —>■ K*'^£~i'e decay can be used to unambiguously verify the left-handed 
nature of the transition operator as predicted by the SM. We provide numerical estimates for the 
relevant angular observables and the resulting decay distributions on the basis of available form- 
factor information from lattice and sum-rule estimates. In addition, we pay particular attention 
to suitable combinations of angular observables in the decays Bs —>■ Kn)£~V£ and B —>■ 

K'k )£'^£~, and find that they provide complementary constraints on the relevant 6 —>■ s short- 
distance coefficients. As a by-product, we perform a SM fit on the basis of selected experimental 
decay rates and hadronic input functions, which results in \Vub\ = (4.07 ± 0.20) ■ 10“®. 


I. INTRODUCTION 

The value of |T4f,| represents one of the least 
well-measured parameters in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix of the standard model (SM). 
Moreover, at present, its inclusive determination from 
B —>■ XuiiTi decays and the extraction from exclusive 
semileptonic or leptonic decay modes lead to somewhat 
different results (see e.g. the review in my Indepen¬ 
dent phenomenological information on 6 —> m transitions 
will clearly help to better understand the origin of these 
discrepancies and the underlying theoretical uncertain¬ 
ties. As the solution to this |I4f,| puzzle might also be 
related to physics beyond the SM, one should also take 
into account possible new physics (NP) effects; see [SHI] 
for recent work in that direction. 

The proliferation of unknown parameters, which 
arises in a model-independent approach with generic 
dimension-6 operators in the effective Hamiltonian, can 
be handled with a sufficient number of independent 
experimental observables in & —>■ u transitions. An 
example is H —>■ TnT)ii'£ where the analysis of the 

secondary p tttt decay introduces a large number of 
angular observables with different sensitivities to the 
individual short-distance coefficients |1|. This is similar 
to what has been extensively used in the analysis of rare 
exclusive b s£~^£~ transitions [5Hini- Because of the 
large hadronic width of the p resonance and the question 
of the S- and P-wave composition of the experimentally 
measured dipion final state, a precision determination of 
\Vub\ from this decay also requires a better theoretical 
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understanding of the B —>■ decay spectrum [nnnj. 

In this article, we focus on the decay Bg —>■ A'*+(—> 
K'k)£~V£, which provides similar insight into the short- 
distance couplings as the decay B p(—)• 7r7r)t'“p£. How¬ 
ever, the width of the K*-meson is sufficiently smaller 
than that of the p resonance, P^* ~ Tp/d ~ 50 MeV. 
Moreover, from studies of the decay B —>• K*J/'ip the 
S-wave background below the K* resonance in B decays 
is constrained to small values, with the S-wave fraction 
Fg < 7% on-resonance [T3]. The decay Bg —)• K'*+(—)• 
K'k)£~V£ thus provides a promising alternative channel 
for a precise determination of \Vub\ in the SM, as has 
already been advocated for in M- For the same rea¬ 
son, it can also be used to constrain NP contributions in 
b ^ u transitions, in particular, as we will show below, 
to exclude possible effects from right-handed currents. 

Another benefit of the decay Bg —> K*~^£~p£ is the 
opportunity to combine it with the rare B —> K*t^£~ 
decay. The secondary decay K* —>■ Ktt is identical in 
both decays, which leads to a one-to-one correspondence 
between angular observables. Hadronic form factors in 
both decays are related by the SU (3) / symmetry of the 
strong interaction, and therefore hadronic uncertainties 
in ratios of angular observables from the two decays are 
expected to be under control^ 


^ The B K*£'^£~ decay amplitude also receives corrections from 
nonfactorizable (i.e. not form-factor like) contributions involving 
hadronic operators in the b s effective Hamiltonian. Semilep¬ 
tonic b u transitions are free of such effects. A comparison of 
the two decays can thus also shed light on the size of nonfactoriz¬ 
able hadronic matrix elements and the validity of the underlying 
theoretical framework. A detailed study along these lines is be¬ 
yond the scope of the present work. 
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Furthermore, these ratios of angular observables are 
sensitive not only to \Vub\-, but also to bilinear combina¬ 
tions of the Wilson coefficients describing semileptonic 
b ^ u and radiative 5 —>■ s transitions in the SM and be¬ 
yond. In light of the present deviations between LHCb 
measurements and the respective SM predictions for a 
few angular observables in the B —>■ K* channel (see e.g. 
[m [n], and also HZl), we will show how this can be 
exploited to obtain complementary information on the 
6 —>■ s Wilson coefficients. 

The outline of the article is as follows. In section El 
we introduce the effective Hamiltonian for semileptonic 
b —>■ uiPi transitions, including NP operators, and de¬ 
fine the angular observables for Bg —t iF*(—>■ K7r)£i'g 
transitions. In the following phenomenological section 


III we identify SM null tests among the angular observ¬ 


ables, and derive expressions in a simplified scenario with 
only right-handed NP contributions. We also define op¬ 
timized observables and highlight the synergies between 
the angular observables in Bg —>■ K7r)i i?g a nd 

B —>■ K*{-^ KTr)£^£~. In the numerical section IV we 
first perform a fit of the Wilson coefficients for (V-A) and 
(V-l-A) currents to experimental data for B —)■ 7r~^£~Pi, 
B~ —)■ t~Dt and B —>■ Xu£~V£ decays. On the basis of 
this fit and theoretical estimates for the relevant form 
factors, we then provide numerical predictions for the 
angular observables and partially integrated branching 
ratios for Bg —>■ A'*+(—>■ K7r)£i'g decays, before we con¬ 
clude in section The helicity basis for the Bg K* 
form factors is defined in appendix where we also infer 
the form factor parameters from light-cone sum rule and 
lattice QCD results. The appendices and are dedi¬ 
cated to details on the determination of the hadronic am¬ 
plitudes and the angular observables of Bg K*^£~i?g 
decays within and beyond the SM, respectively. 


II. DEFINITIONS 

A. Effective Hamiltonian for b u£i/i 

We parametrize possible new physics contributions to 
b —>■ u£ui transitions in a model-independent fashion in 
terms of a low-energy effective Hamiltonian, which can 
be written in the form 

Y.^xOx+ h.c.. (1) 

Here the most general set of dimension-6 operators {Ox} 
is given by 

Ov,^ = [uj^Pib] [hfiPL^i] , 

Os,^ = [mb] [£PlM , ( 2 ) 

Ot = [uaf^’'b] [la^^PLUe] , 

where Pi G {Pl,Pr} are chiral projectors, and we have 
restricted ourselves to (massless) left-handed neutrinos 


and ignored the possibility of lepton-flavor violating cou¬ 
plings. (The generalization to more exotic scenarios with 
light right-handed invisible neutral fermions is straight¬ 
forward, see e.g. [3].) Since in the presence of NP the 
notion of Vub becomes ambiguous, we normalize the op¬ 
erators in eq. Q to an effective parameter Vf}}, which 
can be taken, for instance, as the value of Vub that one 
obtains from a global CKM fit within the SM. If NP ef¬ 
fects are small, one would then have Cv,l — 1 (while 
in the SM Cv,l = 1 and Vub = Vfb’ with all other Wil¬ 
son coefficients vanishing). Comparing with reference [3] , 
where the modifications of left- and right-handed quark 
currents has been parametrized in terms of eL,R together 
with a new mixing-matrix V for right-handed currents, 
our conventions are related via 



£l — Cv,L ~ 1) 


/ Vub 


£R — Cy jzi . 


(3) 


B. Angular distribution in Bg —> K*tvi 


The four-fold differential decay rate for Bg -G K*~^£~Dg 
is defined in terms of the dilepton invariant mass , the 
polar angles 6^ and 6k* in the and K* rest frames, 
respectively, and the azimuthal angle (/> between the pri¬ 
mary and secondary decay planes. 


Stt ^ K*£+De] 

3 dg^ d cos 9^ d cos 9k* d(j) 


J(q^, 9 i, 9 k*, 4 >) ■ 


(4) 


It can be expanded in a basis of trigonometric functions 
of the decay angles. We define 


j{q^,9i,9K*,(t>) = JigSiT? 9 k* + Jiccos^ 9 k* 

+ (J 2 s sin^ 9k* + J 2 c cos^ 9k* ) cos 29g 
+ J 3 sin^ 9k* sin^ 9i cos 2(j) 

+ J 4 sin 29k* sin 29^ cos (f) 

-I- Js sin 29k* sin 0^ cos (j) 

-\- (Jes sin^ 9k* + Jqc cos^ 9k* ) cos 9e 
-I- J 7 sin 29k* sin 9^ sin (j) 

-I- Jg sin 29k* sin 29i sin (j) 

+ Jg sin^ 9k* sin^ 9i sin 2(j ), (5) 

with angular observables Ji(^a) = for i = 1,..., 9 

and a = s, c. By construction, the functional dependence 
of the angular distribution eq. ([^ on the angular observ¬ 
ables is identical to the one for B -G V{^ PiP 2 )£~^£~ 
decays in [5], to which we refer for further details. 

Explicit expressions for the angular observables in 
terms of hadronic form factors and Wilson coefficients 
for b —> u£Di in the general operator basis 0 are derived 
in the appendices. 
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III. PHENOMENOLOGY 

For the remainder of this article we restrict our anal¬ 
ysis to vector-like couplings; i.e., we assume Cg ^ = Cj. — 
^T 5 “ 0 simplicity. This leaves us with only two op¬ 
erators for left- and right-handed 6 —>■ u currents, which 
we refer to as SM-I-SM’. We emphasize that with future 
experimental data one can also test for scalar and ten¬ 
sor currents on the basis of the formulae provided in ap¬ 
pendix 


A. Null tests of the SM 

The twelve angular observables Ji as introduced in eq. 
([^ are not independent. Within the SM, they can be 
expressed in terms of four real-valued quantities: |fVp, 
and the three form factors Fx,||,o- This fact can be used 
to define a series of eight null tests that hold within the 
SM: 


J4 = 6V2 lAp ^ uf F^ll Fo , 


and 


r- , 

J5 = 24y2|7V|2^Re{a2} Fq , 

vr 

J63 = 48|fV|2M|^Re{a2}F^F||, 

is = 12 72 |iV|2 Im {(72} Fj_ Fo , 

V9 

i9 = 24|fV|2M| Im{(72}Fj_F||, 


( 8 ) 


(9) 


together with Jgc = Jy = 0 (all relations valid in the 
SM+SM’ scenario). Here, we introduce a normalization 
factor. 


|iVp^ 


3 • 2i07r3M| 


( 10 ) 


4 .J 2 J 3 + Ji - 4i| = 0 , 

8 J 1 , Ji, - 3Ji - 12 J| = 0 , 

JicJqs — 0 , 

16Jl-36Ji-9Jl = 0, 

Jqc = J 7 = Js = J 9 = 0 ■ (6) 

Deviations from these relations are immediate signs of 
physics beyond the SM. This is in contrast to exclusive 
b —)■ si^i~ decays, where such relations are broken by 
nonfactorizing long-distance contributions. 


B. Angular observables for SM+SM’ 

In the SM+SM’ scenario, we obtain a very simple 
structure of the angular observables, which can be ex¬ 
pressed in terms of hadronic form factors (defined in the 
transversity basis, see appendix 0, and three indepen¬ 
dent combinations of Wilson coefficients, 

= \^V,L =*= ^V,r\'^ 1 

— 2(72 = {Cv,L ~ ^V,R)i^V,L + ^V,r)* ’ 

which depend on the absolute values and ICyi^l 

and the relative phase of the two Wilson coefficients 
(the absolute phase is irrelevant in the angular observ¬ 
ables). Notice that is even under parity transforma¬ 
tions (F -O' F), while (72 is odd. Neglecting the charged- 
lepton mass (which is valid as long as raij +; 1), we 
find 


ks = 3 J 2 s = 9 |7V|2 [( 7 + |Fj.|2 + ( 7 ,- |F|| 

Mi 

Jlc = - J2c = 12 |iV|2 ^ ( 7 - |Fo|7 


and A = A(M^,,q^) denotes the usual kinematic 
Kallen function. The normalization |7Vp is chosen such 


that 

dr 

dq"^ 


E 

A=0.J_,|| 


( 11 ) 


= \NfMl 




9 Ml 

F\k + ^\Fo\- 


where the transversity amplitudes are defined in 
appendix 


Beside the decay rate, one can also define the leptonic 
forward-backward asymmetry HpB via the weighted in¬ 
tegral 

1 /■+^ d^T 

, dcos9,sg„(cose,).jj5.j— 

( 12 ) 

In the SM+SM’ scenario, one finds that Ape takes the 
rather simple form 


ApB 


2Re {( 72 } 1AF|| 


(13) 


Note, that the bilinear (72 is unconstrained by present 
experimental measurements of semileptonic b ^ u tran¬ 
sitions. Therefore a measurement of ApB would provide 
complementary information on the Wilson coefficients. 
In particular, the sign of the forward-backward asymme¬ 
try resolves the present ambiguity between Cy versus 
^v,R^ see section IIVI 

Similarly, the traction of longitudinal K* mesons is 
defined as 


Ft. = 


1 


r-l-l 


dV/dq^ J_T 


dcos 9k* (cos 6k* ) 


d^T 

dq^ dcos 6k* ’ 
(14) 
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where LOF^iz) = (5z^ — l)/2. In the SM+SM’ scenario 
this yields 


Fl = 


'1 1^ Ul 


a+|^^^P+ar(|i^||P + ^|FoP) 


C. Optimized observables in SM+SM’ 


It is now possible to construct particular combinations 
of angular observables where the hadronic form-factor 
dependencies cancel (at least partially), and, as a con¬ 
sequence, these observables are sensitive to the short- 
distance Wilson coefficients, only; or vice-versa. 

We begin with observables where the form-factor de¬ 
pendencies cancel. These can be defined in complete 
analogy to what has been discussed in [3, 


+(5) 

Jj rp 


V 2 J 4 

\l —J2c{‘^J2s — J 3 ) 
h 

2 J2c(2 J 2 S + J^) 

Jqs 

2^(2i2c)2 - {JzY ’ 
2J8 

•y/—2 J2 c(2 J 2 S + J^) 

_~^9_ 

y(2J2c)2 - ( 


(16) 


Within the SM+SM’ scenario, the form factors depen¬ 
dencies cancel exactly at every point in the spectrum. 
However, for integrated angular observables one has to 
take into account the different kinematical prefactors, 
and a residual form-factor dependence will remainj^ In 
the SM+SM’ scenario these optimized observables read 


^ = 1 , 

JFl rp - H rp 

= 4^) 


= 2 


= 2 


Re { 0 - 2 } 


r+^- 


Im{(T2} 




(17) 


We continue with the construction of observables that 
are only sensitive to form-factor ratios. Just as in i? —>■ 


^ We emphasize again that the cancellation of form-factor de¬ 
pendencies holds for the whole spectrum, in contrast to 
B —>■ where it can be spoiled by contributions with 

intermediate photons dissociating into . 


K*£'^£~, we find that the SM+SM’ scenario solely allows 
us to extract one form factor ratio, namely Fq/F\\^ in five 
different ratios of angular observables, 

Fo{q'^) _ _ —J 2 C 

v4^ii(9') "47 ”7174 


(18) 

Inconsistencies among these relations would indicate NP 
beyond the SM+SM’ scenario. 


J 4 


-J'. 


2c 


2 j2 .q ~ Jr 


‘2‘J2fs — Jr 


V2Js 

-Jo 


D. Synergies with B —> 


The decay B K*{-^ KTr)£~^i~ is induced by the 
flavor-changing neutral current (FCNC) transition b —>■ 
sl'^l~. At low hadronic recoil, > 15 GeV^, it is again 
dominated by four-fermion operators which can be ex¬ 
tended to a SM+SM’ scenario. The structure of angu¬ 
lar observables JnW'') in those decays is similar as for 
Bg —>■ K*^£vi. The analogous combinations of Wilson 
coefficients which enter the Jn{q^) now read pf and p 2 - 
(For the explicit definition and a detailed phenomenolog¬ 
ical discussion, we refer to [9].) 

With this we can define a number of useful ratios of 
angular observables Jn{q^) in H —)• K*£~^£~ and Jn{q^) 
in Bg —>■ K*~^£Di, 

R„(9’) S AA , (19) 

Jn[q^) 

for n = Ic, 2c, 4, 5, 6s, 8, 9, as well as 


Ri±{q^) 

R2±{q^) 


2Jis(g^) ± SJsjq'^) 

2Ji,(g2)±3J3(92) ’ 

2J2.(g^)± J3(g^) 

24(g2)±J3(g2) ■ 


( 20 ) 


Within these ratios, the dependence on the hadronic form 
factors can be expected to cancel up to corrections from 
the violation of the SU{3)f symmetry of strong interac¬ 
tions, and from nonfactorizing hadronic matrix elements 
in exclusive b —^ s£'^£~ transitions. In the limit where 
these corrections are neglected, we find 

for n = 1+, 2+ 

for n = 1—, Ic, 2—, 2c 


for n = 4, 5, 6s 
for n = 8, 9 

( 21 ) 

Of particular interest are ratios that are proportional 
to the combination p 2 oc Re {C 7 g(g^) C*o}, where in the 


R„ ~ 


87r^ 


\VtbVt 

I Kb 


* |2 
isj 


Pi 


Pi 


Re 1/32} 
Re {cr 2 
Im { 


P2 


. ImlcTol 
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Decay 

]GeV^] 

Measurement 

Reference 


- 

(1.70 ±0.80 ±0.20) • 10-'‘ 

m 

B~ —> t~Ut 

- 

(1.25 ±0.28 ±0.27) • 10-'‘ 

m 

- 

(1.83j:i];g T 0.24) • 10-'‘ 

m 


- 

± 0.11) • io-'‘ 

m 


[0,2] 

(1.280 ±0.196) • 10-® 



[2,4] 

(1.192 ±0.135) • 10-® 



[4,6] 

(1.446 ±0.108) • 10-® 

m 


[6,8] 

(1.437 ±0.105) • 10-® 


[8,10] 

(1.525 ±0.106) • 10-® 



[10,12] 

(1.490 ±0.111) • 10-® 



[0,2] 

(1.173 ±0.219) • 10-® 



[2,4] 

(1.526 ±0.103) • 10-® 



[4,6] 

(1.213 ±0.105) • 10-® 

m 


[6,8] 

(1.465 ±0.102) • 10-® 


[8,10] 

(1.473 ±0.108) • 10-® 


—>■ TT^fJ.~U-r 

[10,12] 

(1.404 ±0.124) • 10-® 


[0,2] 

(1.225 ±0.182) • 10-® 



[2,4] 

(1.277 ±0.128) ■ 10-® 



[4,6] 

(1.274 ±0.109) • 10-® 

m 


[6,8] 

(1.498 ±0.103) • 10-® 


[8,10] 

(1.405 ±0.115) • 10-® 



[10,12] 

(1.617 ±0.104) ■ 10-® 



[0,2] 

(1.95 ±0.32) • 10-® 



[2,4] 

(1.06 ±0.27) • 10-® 



[4,6] 

(1.51 ±0.28) • 10-® 

m 


[6,8] 

(0.97 ±0.23) • 10-® 


[8,10] 

(0.78 ±0.22) • 10-® 



[10,12] 

(1.59 ±0.28) • 10-® 



TABLE I. Summary of the experimental likelihoods for 
branching fractions of the exclusive b ^ u transitions. We 
assume no correlation among the B~ —>■ t~v't data, and use 
the correlation matrices as given in |22l tables XI and XII], 
[23 table III and IV], [H tables XXIX and XXXII] and [23 
table XVII] for the respective data on —>■ jj.~Vfj, decays. 


SM C-jg{q^) is a linear combination of the Wilson coef¬ 
ficients C 7 ® and in 6 —s transitions (see |5] 

for the explicit definitions). Optimized observables in 
B —^ K*£^£~ only allow to access the ratio |C|®/C;^g|, 
whereas the ratios are sensitive to • Ciq. Measur¬ 
ing the corresponding ratios Jn / Jn thus allows to directly 
access the dependence of and to test the theoret¬ 
ical predictions which are based on an operator product 
expansion in the heavy 6 -quark limit. 


IV. NUMERICAL RESULTS 


In this section we derive numerical results fo r the 
angular observables J„ as introduced in section 


IIB 


Our analysis is carried out within a Bayesian frame¬ 
work, for which we use and extend EOS [20] for all 


numerical evaluations. As prerequisites to our numerical 
study of the angular observables, information on the 
Bg K* form factors, and constraints on the b ^ u 
Wilson coefficients are needed. These will be expressed 
through a-posteriori probability density functions 
(PDFs) labelled P(0FF|theory) and P(0AB|exp. data), 
respectively. We refer to appendix for the precise 
definition of P(0ff I theory). 


A. Determination of Cv,l and Cv,r 


For the following numerical analysis we consider ex¬ 
perimental data on the branching ratios for leptonic 
B~ —> t~Vt and semileptonic B^ —>■ decays as 

summarized in table |Tj together with the averaged value 
for \Vub\ from the inclusive determination, 

1 = (4.41 ± 0.21) • 10-3 P . (22) 

Within the SM-I-SM’ scenario, the additional right- 
handed operator contributes differently to the individual 
decay rates, corresponding to (see e.g. P) 

IK’f P \^V,L ~ ^V,r\ ) 

^ IKf P \Cv,l+Cv,r\^ , (23) 

\viT'-\^^\v:^mM+\cv,Rn ■ 

In order to illustrate the NP reach of our analysis, we fix 
the auxiliary parameter to a value that lies between 
the exclusive and inclusive determinations of j V^b \ within 
the SM, 

IKf I = 3.99 • 10-3 . 


With this we can constrain the absolute values and the 
relative phases of the Wilson coefficients and Cy^, 
where the SM-like solution would correspond to jCy ^ | ^ 


1 and C 


V,R 


0 . 


We construct a likelihood P(data|0AB, Af) from 
(multi)normal distributions as indicated in table |l] and 
eq. (22). Note that we assume that the results for the 
B~ —)■ T~Vr branching ratios [H] and [20] are uncor¬ 
related, since the underlying sets of events use different 
tagging methods for the selection process. The same as¬ 
sumption applies to the results of [21] and m- At this 
time, we only use theoretical input from light-cone sum 
rules for the B ^ tt transition form factors, and therefore 
restrict ourselves to the kinematic range q^ < 12 GeV^. 
For a consistent inclusion of lattice results on the P —)■ tt 
form factor in the high-g^ region (see e.g. [281130] . but 
also note added below), we presently do not have access 
to the necessary correlation information required for our 
statistical procedure. 

Within our analysis, we address the theoretical uncer¬ 
tainties using nuisance parameters for the hadronic ma¬ 
trix elements. These are the P-meson decay constant 
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fg -, and the parameters of the the B tt vector form 
factor its normalization as well as two 

shape parameters see [3T] for their definition. For 
the B-meson decay constant we use a Gaussian prior 
with central value and minimal 68% probability inter¬ 
val fg- = (210 ± 11) MeV, as obtained from a recent 
2-point QCD sum rule at NNLO accuracy [35]. As prior 
for the form factor parameters we use the a-posteriori 
distribution obtained from a recent Bayesian analysis of 
the LCSR prediction at NLO accuracy m- 

In order to assess the physical implications of possible 
deviations from the SM expectations, we compare the 
fit results for three different scenarios. In all cases we 
assume Cv.l to be real-valued (i.e. a possible NP phase 
in the left-handed b —>■ u transition should be associated 
to V^^). As already mentioned, the fit to the considered 
data is only sensitive to the relative phase between the 
Wilson coefficients Cv,l and Cvm, and consequently we 
will always encounter an irreducible degeneracy related 
to Cvx/R — t —Cv,L/R- 


1. First, we consider the scenario “left” that fea¬ 

tures only the left-handed current. In this 
case the number of parameters is five, = 

2. Next, we consider the scenario “real”, in 

which Cy g is present and real-valued. The 

set of AB parameters then reads 9'^^g = 

{Cy^g, BeCyg, /f-(0), bf\ b^\ fg-) . 


3. Last but not least, we also consider the sce¬ 
nario “comp”, which includes a complex-valued 
with the full seven parameters, 0‘^g'^ = 
if'V,L^ bf'^, b2^, Jb-)- 

For all scenarios (M = left, real, comp), we obtain the 
a-posteriori PDF as usual via Bayes’ theorem, 


-P(^AB|data, M) 


P(data|0AB, Ar)Po(^AB, Al) 

P(data, M) ’ ^ ’ 



Significance [cr] 



Quantity 

“left” 

“real” 

“comp” 

d.o.f. 

Reference 

/+(t) 

3.11 

2.36 

2.36 

3 

m 


+0.57 

+0.39 

+0.39 

1 

HU 

B~ —>■ T~l'-r 

+0.64 

+0.34 

+0.34 

1 

m 

+0.99 

+0.75 

+0.75 

1 

m 


-1.84 

-2.35 

-2.35 

1 

m 


0.85 

1.08 

1.08 

6 

HD 

—>■ 

0.87 

0.98 

0.98 

6 

m 

1.70 

1.97 

1.97 

6 

m 


2.53 

2.46 

2.46 

6 

m\ 

B —>■ vi 

+1.67 

+1.45 

+1.45 

1 

HI 


TABLE II. Significances of the measurements at the best-fit 
point closest to the SM point for all three fit scenarios. Notice 
that the pull for the LCSR calculation of the B —>■ tt vector 
form factor /+, marked by a j, does not enter the goodness- 
of-fit calculation. 


1. Scenario “Left” 

Our findings for the scenario “left” can be summarized 
as follows. We find two degenerate best-fit points corre¬ 
sponding to \Cy g\ ~ 1. The best-fit point (with positive 
Cv,l) reads 

= (1.016,0.232, -3.163, -bO.425,0.206). (27) 

We find at this point = 18.54, for 28 degrees of 
freedom (from 29 measurements reduced by 1 fit param¬ 
eter). As a consequence, this represents an excellent fit 
with a p-value of 91%. The significances of the indi¬ 
vidual experimental inputs are collected in table |ll| The 
one-dimensional marginalized posterior is approximately 
Gaussian, and yields 

= 1.02 ± 0.05 at 68% probability. (28) 

Equivalently, this result can be expressed as \Vub\ = 
(4.07 ± 0.20) • 10“^ at 68% probability. 


where 

P(data, M) = / d0AB P(data|dAB, Af)Po(0AB, Af) 

(25) 

is the evidence for the scenario M. The likelihood 
P(data|0AB, at) has already been introduced earlier. In 
all three scenarios, we use for the priors of the Wilson 
coefficients uncorrelated, uniform distributions with the 
support —2 < Cj < -1-2. For model comparisons, we nor¬ 
malize the model priors for the various fits scenarios. The 
corresponding relations read 

Po(comp) : Po(real) : Po(left) = 1 : 4 : 16 (26) 


2. Scenario “Real” 

For the scenario “real”, we find a four-fold ambiguity 
in the data; see figure for an illustration. All local 
modes are degenerate. We calculate the goodness of fit 
in the local mode closest to the SM, 

^AB — 

(1.025, -0.079,0.251, -2.884, -b0.196,0.200), (29) 

and obtain Xreni — 20.47. This fit’s p-value of 81% is very 
good. However, note that the value has increased in 
comparison to the previous scenario. This result war¬ 
rants a comment. The additional degree of freedom in 
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FIG. 1. (left) Contours of the 68% (dark orange area) and 95% (orange area) probability regions for the Wilson coefficients 
Cy and Cy as obtained from our fit. See the text for details. Overlaid are the 68% and 95% contour lines for 
(blue solid lines, negative slope), B~ —>■ £~V(^ (blue solid lines, positive slope) and inclusive B —>■ Xui~Vi (green solid rings). 
The black diamond marks the SM point, (right) Contours of the 68% and 95% probability regions for the Wilson coefficients 
(solid orange lines) overlaying the 68% (dark gray area) and 95% (light gray area) probability regions as obtained from a 
hypothetical measurement of Afb = ± 10 %. 



FIG. 2. Contonrs of the 68% (dark orange area) and 95% (orange area) probability regions for the Wilson coefficients Cy i^ 
and Cyn as obtained from our fit in scenario “comp”. See the text for details. The black diamond marks the SM point. 











































































































form of Cy allows the fit to move the form factor pa¬ 
rameters bi and 62 closer to the central values of 
the prior. This shift occurs at the expense of increas¬ 
ing the significances of the experimental data, while si¬ 
multaneously reducing the significance of the nuisance 
parameters. For completeness, we also list these signifi¬ 
cances for all scenarios in table El The one-dimensional 
marginalized posterior distributions for this scenario are 
approximately Gaussian and symmetric under exchange 
Cyj^ O ReCy ^. We find (at 68% propabality) 

= 1.02±0.05 and | Re < 0.10, (30) 
or 

iReCy^l = 1.02±0.05 and |Cy_i|<0.10. (31) 


with purely left-handed 6 —> u currents over the other 
scenarios very strongly. Moreover, the scenario “real” is 
substantially favored over the scenario “comp”. 

This means that - despite the observed tensions be¬ 
tween the different SM determinations of \Vub\ - a NP 
scenario with right-handed currents does not lead to a 
more efficient description of the experimental data. We 
emphasize again that the statistical treatment of the the¬ 
oretical uncertainties on the hadronic input parameters, 
which are still relatively large at present, has been crucial 
for this argument. On the other hand, the experimental 
data on the inclusive and exclusive decay rates alone also 
cannot exclude large right-handed currents. 


3. Scenario “Comp” 


B. Predictions for Angular Observables J„ 


We repeat the fit in scenario “comp”. As a consequence 
of the additional degree of freedom, the four solutions 
from the previous scenario now become connected. This 
is illustrated in figure We calculate the goodness of fit 
in the local mode closest to the SM, which now reads: 

^conip,!i: _ 

(1.025, -0.080,0.000,0.251, -2.885, -b0.196,0.200). 

(32) 

The individual significances are listed in table |TE and 
amount to a total = 20.48. For the increase of with 
respect to the “left” scenario, see our earlier comment. 
With 26 degrees of freedom the p-value is 77%, which 
is still very good. It is not sensible to provide the 68 % 
probability interval of the one-dimensional marginalized 
posterior, since the solutions are strongly connected. We 
show the contours of the probability regions at 68 % and 
95% probability in figure 

4 . Comparison 

We proceed with a comparison of the various fit sce- 


narios by 
calculated 

means of the posterior odds. The latter can be 
as 


P(Mi data) P(data Mi) Po(-^i) 
P(M 2 data) P(data M 2 ) Po(M 2 ) 

(33) 

We find 

P(“left”|data) . 

P( “real” data) ’ 

(34) 

and 

P(“rear|data) _ ^ 

P(“comp” data) 

(35) 


Using Jeffrey’s scale for the interpretation of the posterior 
odds |33] , we find that the data favour the interpretation 


We can now proceed to produce predictive distribu¬ 
tions for the angular observables in Bg —>■ A'*+(—> 
KTr)£i'i, for which we have two main applications in 
mind. 


1. SM Scenario 


First, we assume the SM case; i.e., we go back to 
—)■ Vub with Cy^ = l and C^ = 0. In this case, 
only the a-posteriori PDF on the Bg —>■ K* form fac¬ 
tors is needed. We obtain the joint posterior-predictive 
distribution for the angular observables by means of 


P{J) = / d^FF I^(fi*FF|theory) J(J — J(0 ff)) ■ (36) 


In practice, the above is carried out by calculating the J„ 
for a set of samples drawn from the a-posteriori PDF. 
In our analysis 10® samples are used. Our results for 
the angular observables, normalized to the decay width, 
are compiled in table m We single out the branching 
ratio, which appears to be the most immediate candidate 
for upcoming measurement. We present our results in 
units of I Tub I which is convenient to extract |Fub| from 
future data. Our results read 

f&GeV^ UR 

/ dg^—= (5.08^°;®®)|Kb|-^ 

GeV2 dg 

dq^ ^ = {8.50f°il)\Vub\-\ 
dg2 ^ = (27.25™) IKbl-^ (37) 

In the above, g%„ = 0.02, and = {Mb^ — 
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2. SM+SM’ Scenario 


Second, we consider the interesting prospect of NP ef¬ 
fects entering the b ^ u transitions, which, according to 
the discussion in the previous subsection, cannot yet be 
ruled out. Based upon our model comparison, we choose 
to give predictions for the scenario “real” only. In order 
to investigate the NP effects on the angular observables 
in Bs —>■ we compute the joint predictive distri¬ 

bution that arises from both posteriors P(^AB |dat a) and 
P(0ff I theory). Our findings are listed in table IV for our 
three nominal choices of bins. In addition, we find for 
the partially integrated branching ratios in the scenario 
“real” 


Jl GeV= 

r^marx. 

/l4.18GeV2 
7o. 02 GeV2 


10 


-5 


^ = (9-5 =b 1.9) 

1 YJ 

^ = (1.55 ±0.19) • 10 


-4 




dq^ 

dB 

dq^ 


= (4.92 ±0.69) • 10 


-4 


(38) 


We also consider suitable ratios of partial decay widths 
over either the —>■ or 


in Bs 
B- - 


T Vr widths. We define three such ratios. 


Rn = 


Ri\ = 


i?, = 


/■ 


dq^lA^l^ 


T{B —>■ T Vr) 
dg2 \A\\^ 

T{B- T-v^) ^ 9r(P- ^ t-D^) ’ 



r^max 


dq^\Al\^ 


(r(po ^ 

8Jis 


I'n)) 

±12J3 


9(r(po 


r+/i i 


4) 


(39) 


where, as already explained above, we only use the 
LCSR-accessible part of the —)■ phase space. 


(r(P° 7r+/r v^)) 


flj 

12GeW dr(P° 


dg2 


TT+fl U^) 


dg2 


(40) 


The ratios Po,||.± ^^e independent of NP effects in this 
scenario. We find numerically, 

i?o = 2.00;°;f2, 

i?ll = 1.36;°;il, 

Rl = 0.79;° j^ , (41) 

where the uncertainties are purely due to the imprecise 
theoretical knowledge of the Bg —t K* form factors, the 
B ^ TT form factors, and the P-meson decay constant. 
Here, correlation information among the various hadronic 
matrix elements would help in reducing these uncertain¬ 
ties. 


n 

(a) 

(^n)/(r) 

(b) 

(c) 

Is 

0.144 

0.368 

f\ OQQ “1" 0.018 
U.Z09 0.020 

Ic 

0.558 

0.260 lSJ:E)?g 

0 Q7Q "1" 0.026 
U.OZO _ 0.024 

2 s 

0 

0 

00 

1 + 

00 

bb 

00 

-41-41 

n 1 OQ b 0.003 
U.IZO _ 0.002 

0.094 lg:gg? 

2 c 

-0.5581 

-0.260 

n Q7Q 0.024 

— U.O/O 0.026 

3 

- 0.010 

-0.129 

0 001 + 

— U.UDi _ 0.009 

4 

0.168 

n 990 0.003 

U.ZZU _ 0.003 

0.19818:gg| 

5 

-0.304 

-0.242 t°o:°o°ol 

-0.29418:[]J§ 

6 s 

-0.189 

-0.4071 

-0.346 lg:gi! 


TABLE 111. Estimates for the normalized nonvanishing an¬ 
gular observables in the SM. The integration ranges are (a) 
IGeV^ < < 6GeV2, (b) 14.18GeV^ < q'^ < 19.71 GeV^, 

and 0.02 GeV^ < < 19.71 GeV^. We normalize the in¬ 

tegrated angular observables (Jn) to the partially integrated 
decay width (T) for the same integration range. 


n 

(a) 

(^n)/(r) 

(b) 

(c) 

Is 

0 1 Q9 "1" 0.025 
U.lOZ _ 0.017 

0.362 l°:°;g 

A 979 + 0.021 

U.Z Z Z _ 0.021 

Ic 

0.574 lg:S!ii 

0.268 l[j:jjjf 

0 

CO 

00 

1 + 

00 

bb 

to to 

0000 

2s 

coco 
00 
00 
bb 
-h 1 

0 

0 

A 1 91 + 0.003 
U.lZl 0.003 

0.091 lg:gg? 

2c 

- 0.574 iK:Ki 

-0.268 ISJigjJ 

1 

0 

CO 

00 

1 + 

00 

bb 

to to 

0000 

3 

- 0.022 lg;gJi 

A 1 ri + 0.018 

—U.iOi 0.016 

-0.08218:[]?i 

4 

0 1 71 "1" 0.009 
U.i / i 0.009 

A 990 + 0.007 
U.ZZ5 0.008 

0.20718:[]gi 

5 

0 971 "1" 0.033 
— U.ZZ 1 _ 0.036 

A 091 b 0.023 
—U.ZZl _ 0.019 

-0.264 lg:!ji 

6 s 

n 1 70 "1" 0.028 

— U.i / Z _ 0.031 

-0.370 l[j:jji 

A 01 0 b 0.031 
— U.OiZ _ 0.041 


TABLE IV. Estimates for the nonvanishing angnlar observ¬ 
ables Jn in the SM-I-SM’ basis for real-valued Wilson co¬ 
efficients. Constraints on the Wilson coefficient are taken 
from data on exclnsive semileptonic fe —> u transitions, see 
text. The integration ranges are (a) 1 GeV^ < q^ < OGeV^, 
(b) 14.18GeV^ < q'^ < 19.71 GeV^, and 0.02GeV^ < q'^ < 
19.71 GeV^. We normalize the angular observables to the par¬ 
tially integrated decays width (T). Note that the quoted sign 
for the angnlar observables J 5 and Jes corresponds to the SM- 
like solution ( |30[ ) with dominating left-handed current. For 
the solution ( |31| ), one simply has to flip the sign of J 5 and 

J 6s ■ 


V. CONCLUSIONS 


The angular analysis of exclusive Bg —l A'*+(—l 
KTr)£i'i decays provides a powerful tool to measure the 
Cabibbo-Kobayashi-Maskawa (CKM) element |14b| in 
the Standard Model (SM) and to constrain new physics 
(NP) contributions to the underlying semileptonic b —>■ 
uiPi transition. In this article, we have identified rela¬ 
tions among the angular observables that serve as null 
tests of the SM. Furthermore, we have constructed op¬ 
timized observables where, also in the presence of NP, 
the dependence on either the hadronic form-factor or the 
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short-distance coefhdents drops out. The fact that the 
same secondary decay, K* —)• K-k, is used for the an¬ 
gular analysis of the rare B —>• K*i~^£~ decay can be 
phenomenologically exploited by measuring certain ra¬ 
tios Rn of angular observables from both decays. In 
the limit where nonfactorizable effects in B ^ 
as well as SU{Z)f symmetry corrections to form-factor 
ratios can be neglected, the ratios i?„ are only sensi¬ 
tive to short-distance coefficients. In particular, we have 
shown that in this way one can directly access the q^- 
dependence of the effective Wilson coefficient function 
in B —>• transitions. 

We have combined presently available experimental 
data on inclusive and exclusive, leptonic and semilep- 
tonic b ^ u transitions with theoretical information on 
hadronic form factors and decay constants, thereby ob¬ 
taining detailed numerical estimates for angular observ¬ 
ables and partially integrated decay widths in Bg —)■ 
Ktt)£D£. Here, we also allowed for the presence 
of right-handed currents that could arise from physics 
beyond the SM. Using a Bayesian approach for the sta¬ 
tistical treatment of theoretical uncertainties, we have 
found that - despite the present tensions between dif¬ 
ferent \Vub\ determinations - the SM is still more effi¬ 
cient in describing the experimental data than its right- 
handed extension. In a simultaneous SM fit to —>■ 

(using light-cone sum rule results for low dilep¬ 
ton mass), B~ —>■ t~Dt and B —>■ X^lvi data, we find 
\Vub\ = (4.07 ± 0.20) • 10“^ with a p-value of 91%. 

On the other hand, right-handed contributions cannot 
be excluded, either. In a SM-like scenario with domi¬ 
nating left-handed currents, we found that the ratio of 
right-handed over left-handed currets is constrained to 
< 10%. Since the decay rates alone are invariant under 
parity transformations, a second solution, with the role 
of left- and right-handed quark currents interchanged, 
is always presents Again, some of the angular observ¬ 
ables in Bg —)■ Ar*^(— )■ KTr)i9£, e.g. the leptonic forward- 
backward asymmetry, are “parity”-odd and can thus un- 
ambigously test the (dominating) left-handed nature of 
semileptonic 6 —>■ u currents. In this case, one would 
obtain strong constraints on the flavour sector of NP 
models with generic right-handed currents. (For a recent 
attempt to construct a left-right symmetric NP model 
based on the Pati-Salam gauge group, which can acco¬ 
modate naturally small right-handed b ^ u currents, see 

EH) 

A crucial ingredient of our analysis has been the im¬ 
plementation of hadronic uncertainties. Improvements of 
our theoretical understanding of nonperturbative QCD 
effects (see also notes added below) would lead to more 
stringent constraints on the value of |t46| and the possible 
size of right-handed b ^ u currents. In particular, predic¬ 
tions from lattice or light-cone sum rules for form-factor 


® Notice that the lepton current with a light SM-like neutrino is 
always considered to be left-handed, only. 


ratios with B and Bg initial states (including correla¬ 
tions between input parameters), and similarly between 
B —TT form factors and the B-meson decay constant, 
would be helpful in this respect. 


Notes added: 

In the final phase of this work, the LHCb collaboration 
measured the ratio of the exclusive semileptonic branch¬ 
ing fractions of Af, —pfjT and A;, —Ac^“ i>^ [351ISS] ■ 
Assuming SM-like b —>■ transitions, with knowl¬ 

edge of the magnitude of \Vcb\ and using information on 
the relevant form factors [37], this ratio can be used to 
extract the branching fraction B{Kb —)■ ppL~v^). As such, 
the branching fraction is a very powerful new constraint. 
However, in light of the present tension in the determi¬ 
nation of Vcb from both inclusive and exclusive b —>■ c£vi 
decays, and in order to follow the logical line of this arti¬ 
cle, the new LHCb measurement should only be used in 
a setup that accounts for NP in both b ^ u and b ^ c 
transitions. 

Another article |5B] that was published in the final 
phase of this work provides updated LCSR results for 
the hadronic form factors for Bg —>■ K* transitions, which 
include correlation information among the form factors. 
This development will help to further reduce theory un¬ 
certainties for this decay. 

In recent lattice studies of the B —>■ tt form factors 
[39l [40] , also the correlation matrix between the relevant 
hadronic fit parameters has been provided. This will also 
allow to include the high-q^ data for the B —>■ decay 

in our statistical procedure, which could and should be 
used in future updates of our results. 
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Appendix A: Form Factors 

There are in general 7 independent hadronic form fac¬ 
tors for Bg —>■ K* transitions. Commonly, these are de¬ 
noted as U, Aq i^ 2 ^ 7i. 2.3; see e.g. the definition in m- 
For our purpose, it is more convenient to start with a 
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definition of form factors in a helicity basis, 

Po= -^^-{K*{k,r])\uf^{l - j5)b\Bs{p)), (Al) 
Pt= ^|-(is:*(fc,? 7 )|M^*(l- 75 ) 6 |i?s(p)), 
and 


q^ [GeV^] 

0 15.00 19.21 

V{q^) 

Ai(g^) 

A2(g") 

AMq^) 

0.311 ±0.026 0.872 ±0.066 1.722 ±0.062 

0.233 ±0.023 0.427 ±0.015 0.548 ±0.015 

0.181 ±0.025 

0.342 ±0.016 0.408 ±0.016 


q^ [GeV^j 

V 

15.00 19.21 

Al 

15.00 19.21 

Ai2 

15.00 19.21 

15.00 

19.21 

1.000 0.271 

- 1.000 

1.000 0.305 

- 1.000 

1.000 0.334 

- 1.000 


P± = {k,r])\uaf,^e'^q''{l + -f5)b\B,{p)), 

^ f^ {K*{k,p)\uafj,^e^*q''{l + -/5)b\B,{p)) , 

Ms^Vq 

(A2) 

which is related to the one proposed in (42]. However, 
compared to |42) . we have chosen a normalization con¬ 
vention such that all form factors are finite in the limit 
—>■ t_ = {Mb, — Mk*Y, and nonzero in the limit 
q^ —)• 0. In the above definition, p denotes the physi¬ 
cal polarization of the K* meson, and e stands for an 
auxiliary polarization vector of the dilepton system with 
polarization states t,±l,0. Notice, that the form factor 
for the pseudoscalar current is not independent, but from 
the equations of motion can be related to Ft, 

Ml 

{K*{k,p)\u'y5b\Bs) = -i -— Ft. (A3) 

mb -f mu 


Instead of the helicity form factors F±, we will use the 
linear combinations 

(A4) 

which simplify the analytical expressions for the angular 
observables. The explicit relations between our and the 
traditional form factor basis read 




Mb, {Mb, + Mk* ) 




for the vector form factor. 


(A5) 


F. = ^ A ,, 


Fn = 


Mb, 

{Mb, + Mk^ )2 {Ml - Ml, - ql Al - A A2 


2Mk*MI {Mb,+Mk>) 


Ft = 


8Mk* Ai 2 
Mb, 

^ A 


Ml, 


(A 6 ) 


TABLE V. Theory inputs for the Bg —>■ K* form factor fits. 
Top: Form factor values at = 0 are taken from LCSR cal¬ 
culations in ^T]; values at — 15 GeV^ and = 19.21 GeV^ 
are taken from lattice QCD simulations [43|. Bottom: Gor- 
relation information for the lattice QGD inputs. The lattice 
QCD values and correlations are produced from the joint PDF 
given in [43 [ Table XXIX] using 5-10® samples. 


for the axialvector currents, and 


pT _ rp 


j^T ^Ml^ Ml.) ^ 

^11 Ml^ ’ 

{Ml^-Ml.){Ml^+3Ml.-qlT2-XT3 

° 2Mk> Mb, {Ml^ - Ml.) 

_ AMk* T 23 

Mb, + Mk* ’ 

(A7) 

for the tensor current. In the above equations, the form 
factors Ai 2 and T 23 are defined as in [43]. 

The form factors fulfill endpoint relations |42l I44j^ 
which in our convention read 


lim Fj_ = lim Ft = 0 , 

F|| ^Mb, 

lim — = —--;-, 

Fq Mb, — Mj^r 

with t± = {Mb, ± Mk^Y- We will use these relations 
for our form-factor parametrization in the numerical fit. 
To this end, we consider a modified “z-expansion” and 
write 


X [ 1 -h5j_ (z(g^^o) - 2(0,<o))] , 

F||.o(9")=F(g^M|JF||,o(0) 


Note that the endpoint relation for the ± form factor in |42[ 
appendix B] should read limg 2 _^^ _ Bv,\IBv,2 = 0. 
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X [l + 6||,o (z(g^^o) - 2:(0,to))] , 

X [l + bt{z{q^,to) - z{0,to))] . (A9) 


e.g. gSMl) 


, u\ - - a - \/4 - b 

zia, b) = ---■ ^ - _. 

y/t^ — a + y/tj^ — b 


(AlO) 


Here, the prefactors contain global kinematic factors, the 
form-factor normalization at = 0 , together with the 
leading pole behaviour from the lowest resonances above 
the semileptonic decay region, P[q^,M‘^)~^ = l — q^/M"^. 
The remaining g^-dependence for each form factors is 
parametrized by a shape parameter bi. The variable 
z(g^,to) is obtained from the conformal mapping, (see 


Here we choose Iq = — y/t^{t^ — t_) which minimizes 

\z\ in the decay region. For the resonance masses we 
use Mb = 5279 MeV, Mb* = 5325 MeV and Mb, = 
5724 MeV [48]. The above parametrization eg. ( A9) au¬ 
tomatically fulfills the end-point relation eg. (A 8 ) for F±. 
The end-point relation for Fii/Fq is fulfilled by imposing 


z( 0 ,to) - z{t-,to) 


(, F\m 

\ Fo{0) 


We fit the Bs —>■ K* helicity form factors Tj_,||,o to 
the nine constraints listed in table IS Our fit uses five 
parameters. 


Off = (F_L(0),F||(0),Fo(0),5i,&||) (A12) 


which represent the three normalizations Tj_ y 0 ( 9 ^ = 0 ), 
as well as two independent shape parameters 6_l,||. As 


a-priori probability Po{Off) we choose uncorrelated uni¬ 
form distributions with a generous support (to be com¬ 
pared with (A15) below). 


0<F±\\ o(0) < 1, -10 < 6 _l < 0, -5 < 6 || < -k5 . 

(A13) 

The likelihood P (theory I^ff) is constructed as the 
product of uncorrelated Gaussian likelihoods for each of 
the LCSR results for the form factors V, Ai and A 2 , 
as well as the joint multivariate Gaussian likelihood for 
the lattice QGD results. All of these are listed in table |V| 


The a-posteriori PDF is obtained as usual via Bayes’ 
theorem. 


P(0ff I theory) 


P (theory I 0 ff ) Pq (^FF) 

/ d^FF P(theory|0FF)Po(^FF) 


(A14) 


For all applications here and in section [ivl we draw 10® 
samples from the a-posteriori distribution. 

The best-fit point, and the ID-marginalized minimal 
intervals at 68 % probability are found to be 


Pj.(0) = 0.349 ±0.037, &j_ =-4.9^1;? , 

P||(0) = 0.379 ±0.031, 6|| = ±0.07 ±0.40. (A15) 

Po(0) = 0.314 ±0.041, 


Although the ID-marginalized distributions are symmet¬ 
ric and resemble Gaussian distributions, we find that the 


I t_ 

2 M 2 _ 


[1 ± ^||(^(t-,lo) - 2:(0,to))] 


(All) 


distribution in eg. (A14| is distinctly non-Gaussian. We 


therefore use the posterior samples to carry out the un¬ 
certainty propagation. 


Appendix B: Ba K*{-^ Kf) 1 vi Decay Amplitude 


In this appendix we give details on the parametrization 
of the matrix element for the decay Bg —>■ K*^£~D(, with 
the subseguent decay K*~^ —)• (Ar 7 r)+. We decompose the 
matrix element as in 


M=F{Xs[ir'] +Xp[i^5i^] 

± [hr.l5v] + Xl^f } (Bl) 

with the prefactor 

F = iV2GFVub gK*KTv Dk* \kKF\, (B2) 


and I^rfI = y/ A(M|-,, M|., M^)/2Mk* ■ In the small- 
width approximation we replace the K* resonance by 


\DK*{k^W 


{k^ - MI,Y + MlSl, 


Mk*Tk- 


-5{k^-Ml,) 


(B3) 


where Fr:» denotes the total decay width of the K* me¬ 
son. Since Fr-. ~ r[Ar* — K-k] to very good approxima¬ 
tion, we use 


F 


K* 


|giC*-»R'7i-P|fcRFp 


(B4) 


Our parametrization of the hadronic matrix element 
of P —>■ V{—>- PiP 2 )£~i'£ decays differes from the one in 
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[S] due to different conventions for the Levi-Civita tensor, 
the phase convention for the polarization vectors, and the 
fact that in this decay only left-handed lepton currents 
contribute. We use 


NXs = ^ cos = -NXp , (B5) 

NX^ = -NX^ 

= -1-^ cos 9ve^{0)AQ 
+ ^-smeve^i+)e+^'^[Ai+A^] 

+ ^sm0ve>^{-)e~'-‘^[A^-A^] , (B6) 

(B7) 

NX!^'' = cos6lye^(-h)£:'^(-)^||_L 

- !^£/^(0)e^(-)e-^'^Ao|| . (B8) 


and 


A||i = + 81 VMB,C^Fo^(g 2 ), 


^t± — 

^o|| = 


/\/f2 

4V2N^CrF]:iq^), 

/Vf2 

4V2N^CrF{iq^). 


(BIO) 


Appendix C: Angular Observables for B Vive 
In the limit me —>■ 0, the angular observables Jn read 

>715 = ^ [3|^il^ + 3|A|^|^ + 16|Ao|| 1^ -I- 16|At_L|^] , 

Jic=l[\A^\^ + 2\Aff + 8\A\\pf], 

•hs = + l^f r ~ 16 |Ao|| P — 16 |At_Lp] , 

[\Ao\'^ - 8|^||j.P] , 

•^3 “ g [l^iP ~ + 16|^o||P ~ 16|At_Lp] , 

(Cl) 


Using the normalisation constant N as given in eq. (10) 
and the general operator basis ([^, we obtain for the in¬ 
dividual amplitude contributions 


= -m 




i^V,L ~ Cv,r) Foig"^) , 


Ai = +4 AMb, {Cyp + Cyp) Fp{q ^), 

A\ = -ANMb, {Cyp - Cyp) U||(g2) , 


Af = -4A 

Ml 


meMs^ 


i^v,L ~ ^v,r) 


"I ~{^s,L ~ ^s,r) , 


mb 


(B9) 


■/s = Re {A^Ai + 2y2Ao||Af*} , 

J6« = ^Re{AfAi*} , 

J6c = -6Re{A||^Af*} , 

Jt = ^ Im [a^AI* - 2V2AepAl*} , 

•^8 = ^Im{A^Ai*}, 

J9 = ^Im{AiAf} . (C2) 
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